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1. INTRODUCTION 


The effects of atmospheric turbulence in both horizontal and 
near-horizontal flight, during the return of the Space Shuttle, are impor- 
tant for determining design, control, and "pilot- in-the- loop" effects. A 
non-recursive model (based on von Karman spectra) for atmospheric turbu- 
lence along the flight path of the Shuttle Orbiter has been developed which 
provides for simulation of instantaneous vertical and horizontal gusts at 
the vehicle center-of-gravity, and also for simulation of instantaneous 
gusts gradients. Based on this model the time series for both gusts and 
gust gradients have been generated and stored on a series of magnetic 
tapes which are entitled Shuttle Stmulation Turbulence Tapes (SSTT). The 
time series are designed to represent atmospheric turbulence from ground 
level to an altitude of 120,000 meters, 

A description of the turbulence generation procedure is provided 
in Secti'on 2, The results of validating the simulated turbulence are 
described in Section 3, Conclusions and recommendations are presented in 
Section 4 with Section 5 containing references ctted. Appendix A contains 
the tabulated one-dimensional von Karman spectra while Appendix B provides 
a discussion of the minimum frequency stmulated. Appendices C and D present 
the results of spectral and statistical analyses of the SSTT. A more detailed 
description of the proper use of the tapes is provided elsewhere [1]. 



2. TUFBJLENCE GENERATION PROCEDURE 


The non-recursive turbulence model used to generate the SSTT i s 
based on von Karman spectra with finite upper limits corresponding to the 
dimensions of the Space Shuttle, relative to the scale of turbulence in the 
atmosphere. Because the scale of turbulence increases with altitude while 
the dimensions of the Space Shuttle are fixed, the finite upper limits of 
the von Karman spectra increase with altitude. In order to take into 
account the resulting spectral changes, the atmosphere, extending from 
ground level to 120,000 meters, was divided into six altitude bands. The 
subsections which follow provide a description of the development and 
application of the turbulence generation procedures. 

2.1 BACKGROUND 

The current turbulence model represents the results of the 
development and evaluation of several different turbulence simulation 
techniques. Two of the earlier techniques, which were given serious 
consideration, warrant further discussion. 

Initial efforts involved refinement and evaluation of a turbulence 
model, TBMOD [2], which had been developed elsewhere [3]. This model was 
based on discretization of the Fourier integral representation of turbulence 
and was designed for use with von Karman spectra. Several problems were 
encountered with TBMOD, both theoretical and practical. First, from a 
theoretical standpoint, the assumptions used in the development of the 
shear (gust gradient) simulation were difficult to justify. Second, from 
a practical standpoint, the output of TBMOD, representing turbulent gusts, 
when subjected to Fast Fourier Transform (FFT) spectral analysis, did not 
possess the proper von Karman spectral shape. Because of such problems, 
further development of the model was halted. 

The second simulation technique was based on digital filter 
theory coupled with a combined von Karman - Saffman spectral model [4]. 
Meromorphic functions were used to approximate the various spectra. Based 
on z-transform theory, recursive difference equations were derived from 
such approximations. Such difference equations were then used to generate 
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the appropriate turbulence gusts and gust gradients. Unfortunately the 
recursive difference equations, i n addition to being somewhat complex, 
proved numerically unstable and could not be used i n their original 
recursive form [5] . 

A non-recursive version of the same difference equations was 
subsequently developed in an effort to overcome the stability problem [5]. 
These difference equations were also quite complex but resembled i n form 
the output of digital filters, characterized by some impulse response 
function, with a white noise input. The present model, as described in 
subsections 2.2 through 2.4, was an outgrowth of this resemblance. 

2.2 SELECTION OF ATMOSPHERIC BANDS 

The standard deviations (a^, a 3 ) and the integral scale 
lengths (Lj, Lg, Lg) of atmospheric turbulence are functions of altitude, 
as shown in Table 2-1. Notice should be taken that the values’for cr.. and 
Lj presented i n this table are consistent with those presented i n JSC 7700 
[ 6 ]. Based on the variation of ck and presented in Table 2-1, the 
atmosphere was divided into six altitude bands as presented i n Table 2-2. 
Within each band, as also indicated in Table 2-2, characteristic integral 
scales of turbulence were selected for use i n calculating the finite 
upper limit of the turbulence spectral model discussed in subsection 2.3. 

2.3 DEVELOPMENT OF VON KARMAN SPECTRA WITH FINITE UPPER LIMITS 

As developed previously [5] the basic three-dimensional von Karman 
relation to be integrated for the dimensionless gust spectra is, 


1 ^1 ’^2 ’^ 3 ^ = 


55 

36air* (l+n 2 ) 17/6 


( 2 - 1 ) 


The corresponding von Karman relation for dimensionless gust gradient 
spectra i s 

s 2(a2-n.2) 

*H/Jj(ra>n2’G3> . 35X3 (m 2,17/6 ^ 

These three-dimensional spectral relations must be integrated over certain 
ranges of values of fig ar| d ^3 to obtain one-dimensional spectral models 

and 
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TABLE 2-1. VARIATION OF STANDARD DEVIATION 
AND LENGTH SCALE WITH ALTITUDE 


ALTITUDE 

(m) 


STANDARD DEVIATION 
OF TURBULENCE 

INTEGRAL SCALES 
OF TURBULENCE 

as 

^(m/sec) 

a 3 (m/sec) 

LjCm) 

L 2 (m) 

L 3 (m) 


1.67 

1.15 

21 

11 

5 


1.98 

1.46 

33 

19 

11 

2.75 

2.20 

1.71 

43 

28 

17 


2.36 

1.89 

52 

35 

23 

2.98 

2.49 

2.05 

61 

42 

29 

3.07 

2.61 

2.19 

68 

49 

35 

3.15 

2.71 

2.32 

75 

56 

41 

3.22 

2.81 

2.43 

82 

63 

47 

3.28 

2.89 

2.54 

89 

69 

53 

3.33 

2.97 

2.64 

95 

75 

59 

3.72 

3.53 

3.38 

149 

134 

123 

31.95/4.37 

3.95/4.37 

3. 9514.3 9 

1961300 

190/300 

1921300 

4.39 

4.39 

4.39 

300 

300 

300 

4.39 

4.39 

4.39 

300 

300 

300 

4.39 

4.39 

4.39 

300 

300 

300 

4.39 

4.39 

4.39 

300 

300 

300 

4.39/5.70 

4.39/5.70 

4.39/5.70 

300/533 

300/533 

300/533 

5.70 

5.70 

5.70 

533 

533 

533 

5.70 

5.70 

5.70 

533 

533 

533 

5.70/5.79 

5.70/5.79 

5.7015.7 9 

533 

533 

533 

5.79 

5.79 

5.79 

533 

533 

533 

51.79/5.52 

5.79/5.52 

51.79/5.52 

533 

533 

533 

5.52 

5.52 

5.52 

533 

533 

533 

5.52 

5.52 

5.52 

533 

533 

533 

51. 52/5.27 

5.52/5.27 

5.52/5.27 

533 

533 

533 

5.27 

5.27 

5.27 

533 

533 

533 

5.27 

5.27 

5.27 

533 

533 

533 

5.27/4.22 

5.27/4.22 

5.27/4.22 

533 

533 

533 

4.22 

4.22 

4.22 

533 

533 

533 

6.01 

6.01 

4.22 

6691 

6691 

955 


Double entries for a tabulated altitude indicate a step 
change in standard deviation or integral scale at that altitude. 
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TABLE 2-1. VARIATION OF STANDARD DEVIATION 

AND LENGTH SCALE WITH ALTITUDE (Continued) 


RLTTTUDE 

(m) 

STANDARD DEVIATION 
OF TURBULENCE 

INTEGRAL SCALES 
OF TURBULENCE 



tr 3(m/sec) 

L d m > 


L 3 (m) 

27000 

7.00 

7.00 

4.22 

20000 

20000 

1230 

30000 

8.23 

8.23 

4.66 

23533 

23533 

1443 

40000 

12.82 

12.82 

6.09 

36693 

36693 

2231 

50000 

18.08 

18.08 

7.51 

51786 

51786 

3128 

60000 

23.94 

23.94 

8.90 

68623 

68623 

4124 

70000 

30.36 

30.36 

10.28 

87063 

87063 

5208 

80000 

37.29 

37.29 

11.65 

106998 

106998 

6376 

90000 

44.70 

44.70 

13.01 

128338 

128338 

7622 

100000 

52.58 

52.58 

14.35 

151010 

151010 

8941 

110000 

60.89 

60.89 

15.69 

174950 

174950 

10330 

120000 

69.62 

69.62 

17.02 

200000 

200000 

11800 


2.3.1 Upper Limits of Integration 

The upper limits of integration for j = 2 and 3 are calculated 
according to the relation [7] 

= aL./Aj (j = 2,3) (2-3) 

= 1,339 

= integral scale of turbulence associated with the 
spectrum 

= characteristic length of Space Shuttle in the 
« 7 'th direction 

Values of L^ for the six bands are given i n Table 2-2 while the character - 
fstic lengths, 7 .$ are presented in Table 2-3. 


^ijmax 


where 
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TABLE 2-2. SUMMARY OF TURBULENCE PARAMETERS 
IN DISCRETE ALTITUDE BANDS 
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TABLE 2-3. 


CHARACTERISTIC DIMENSIONS 
OF THE SPACE SHUTTLE [8] 


Characteristic 

Length 


Z 


1 


Z 


3 


Magnitude 

(ft) (m]_ Explanation 


39.56 

12.06 

mean aerodynamic chord 

39.05 

11.9 

1/2 wingspan 

10.95 

3.34 

1/2 fuselage thickness 


In the case of special consideration must be given to the 

dimensional frequencies corresponding to the dimensionless limits. The dimen- 
sional frequency limit satisfies the relation 

flmax ~ fi i lmax V/ (2TraL.. ) (2-4) 

where 

V = vehicle velocity 

The maximum dimensional frequency which the Space Shuttle simulators are 
capable of handling i s 4 hertz, Thus any higher frequencies should be ex- 
cluded from the simulation. For this reason the dimensionless frequency 
limit ^-}i max roust satisfy the relation 

®ilmax = min(aL.j/Jlp 27raL^ W v > < 2 - 5 > 

where 

f lmax - 4hertz 


Values of G Umax 


based on Eq (2-5) are included in Table 2-2. 


2.3.2 One-Dimensional Spectra 

There are six spectra of primary interest for turbulence simu- 
lation, as indicated in Table 2-4. Based on second-order numerical inte- 
gration, the six corresponding three-dimensional gust and gust gradient 
spectral relations, as given by Eqs (2-1) and (2-2), were integrated over 
Qg and £1^ (with the appropriate upper limits). The resulting one-dimensional 
spectra for all altitude bands are presented in Appendix A. These spectra 
were used in establishing the inpulse response functions associated with 
digital filter simulation processes described in subsection 2.4. 
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TABLE 2-4. TYPES OF SIMULATED TURBULENCE 


Type 

Corresponding 

Spectrum 

Comments 


*11 

longitudinal gust 

u 2 

CM 

CM 

transverse gust 

u 3 

*33 

vertical gust 


*22/11 

yaw 

Sllg/SXj 

*33/11 

pitch 

3u 3 /3x 2 

*33/22 

roll 


2.3.3 Dimensionless Energy Content 

The total dimensionless energy content of each one-dimensional 
spectra tn each altitude band was established by integrating the corres- 
ponding spectra over the appropriate finite limits , indicated in Table 2-2. 
The resulting energy content Is presented in Table 2-5. As might be expected 

TABLE 2-5. DIMENSIONLESS ENERGY CONTENT 
FOR GUSTS AND GUST GRADIENTS 


ALTITUDE 

SPECTRUM 

BAND 

*11 

*22 

*33 

*22/11 

*33/11 

*33/22 

1 

.6225 

.5010 

.2752 

.5877 

.1525 

.1557 

2 

.8595 

.8560 

.8383 

13.147 

12.171 

12.308 

3 

.8956 

.8952 

.8809 

24.767 

22.643 

22.890 

4 

.9298 

.9296 

.9197 

54.123 

49.527 

50.060 

5 

.9977 

.9953 

.9251 

1740. 

41.71 

95.62 

6 

1.000 

.9973 

.9363 

2309. 

52.08 

391.6 
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the total dimensionless energy content of eacto of the turbulent gust series 
is less than unity. The dimensionless energy content for each gust gradient, 
however, is not limited i n such a manner and range as high has 391.6. For 
both gusts and gust gradients the total energy content increases with altitude 
because of similar increases in the limits of integration, 


2.4 DIGITAL FILTER SIMULATION 

As suggested in subsection 2.1, simulated turbulence, Y(t), can be 
interpreted as the response or output of a control system [9J with double- 
sided response functions, h(t), subject to an input consisting of Gaussian 
white noise I(t). This response can be represented by the convolution 
integral 


Y(t) 



h(x) I(t-r) dx 


(2-6) 


Based on filter theory the double-sided spectrum, ^(si), of the simulated 
turbulence satisfies the relation 

$ dy (%) = H(n 1 )H*(n 1 )4 0I (n 1 ) (2-7) 

where H(Q^) = F[h(t)] 

$ D j(Pj) = double-sided power spectrum for white noise 


Generally the standard deviation of any white noise signal has a value of 

unity. Furthermore in most practical situations the white noise is defined* 

to occur over some interval extending from -£L to +Q, , For this case , 

M lmax lmax ’ 

1 

‘l' = " 




2fi 


lmax 


_ T 
" 2rr 


( 2 - 8 ) 


Actually the term "energy" is not precise when dealing with gust gradients. 

The subscript, i, normally applied to the variables and T, has been 

suppressed in Eqs (2-8) through (2-13) for simplicity. 1 
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where 


T 


By substitution, 


time interval associated with generation 
process (-*/0 |max > 


%(%) a H(!]j)H (Pj) 2tt 


If is limited to real values, 

= ^ (^j) 2 tF 


Rearrangement o f Eq (2-10) yields 

H(^) = Vx 


(2-9) 


( 2 - 10 ) 


( 2 - 11 ) 


Then based on the definition of the inverse Fourier transform, the double- 
sided impulse response function h(t) can be expressed as 


where 


h(t) 


00 

*/Vt 

-00 

oo 

3 zj 008 (ty*) 

<4 


■ 2 


= 2 


2ir 


7 T 


cos(n^t) dfij 



COS(Qjt) dfij 


<()y(fij) = single-sided spectrum of Y(t) 


( 2 - 12 ) 
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The single-sided spectra tabulated in Appendix A correspond to $>y(ftj) 

The discrete version of the convolution'integral given in Eq (2-6) yields 


where 

Y(k) 

+N 

= Z h(j) I(k-j) T 
j“-N 



Y(k) 

= discrete sampled turbulence output 


h(j) 

= discrete double-sided impulse 
function h(jT) 

response 


I(k) 

= discrete sampled white noise i 

nput 


Eq (2-13) represents*the basic, non-recursive relation for the generation of 
simulated turbulence . The impulse response functions h(t) were evaluated 
by means of second-order numerical integration of Eq (2-12) using the six 
spectra from Appendix A. I n carrying out this evaluation some maximum value of 
t must be established, corresponding to the value of N for Eq (2-13), and re- 
placing the infinite limit of Eq (2-6). As discussed i n Appendix B this maximum 
time limit determines the minimum frequency, i max * f° r which the corresponding 
spectrum is accurately simulated. 

The values of dimensionless time increment, T., used for the six 

altitude bands are included in Table 2-2 and are based on the values of ft. „ 

imax 

shown in the same table. Thus the Nyquist generation frequencies fi^g for the 
simulated turbulence correspond to the upper frequency limits for ft^ as com- 
puted by Eq (2-5) for each altitude band, 


* 

In certain references 19,111 to correct for the "effect of digitizing" the 
series represented by Eq (2-13) has been divided by /T. This process can be 
seen to be dimensionally incorrect and actually results from the use of a white 
noise spectrum with unit strength instead of a strength of T/2ir. 
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2.5 EFFECTS OF DIGITIZATION 

The effects of digitization in turbulence simulation have been 
considered by a number of investtgators [9-12], As a result of these studies 
two basic digitization effects have been generally identified. 

The first effect results from the assumption of a white noise 
specrrum with u n it strength instead of unit power, noted in subsection 2.4. 

To correct for such an "effect" the proposed procedure i s to divide the 
series approximation of the convolution integral by vT. This "effect" 
disappears when the white noise spectrum has unit power. 

The second effect involves the tapering of the spectrum of simulated 
white noise, , in the vicinity of the Nyquist generation frequency, 

Seme investigators [10,12] have considered it necessary, because of the taper- 
ing effect, to generate the simulated turbulence time series at a rate from 
four to ten times the rate at which the series will be sampled. 

The second effect arises from the discrete processes associated 
with both the generation and sampling of the simulated turbulence. In the 
case of discrete white noise with unit variance, the time series involved i s 
basically a train of step functions as shown in Figure 2-1. The autocorre- 
lation function of the train of step functions depicted in Figure 2-1 can 
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be shown to be 



(M 1 t g ) 
(M > t q) 


(2-14) 


The corresponding double-sided power spectrum by definition is 


*' DI (n) s fer di (t)i 

Tq sin 2 (flT Q /2) 

** (HT g /2) 2 

x sin 2 (W2n NG ) 
^Vig (OTT/2n NG ) 2 


(2 - 15 ) 


The single-sided version of this power spectrum is shown in Figure 2-2, 



Figure 2-2. White Noise Spectra 



Theoretical white noise, by definition, is characterized by a 
uniform power spectral 'distribution. In order to avoid infinite power, such 
a spectral distribution is normally restricted to the frequency band 

c 1 ^ng) f° r a double-sided spectrum. For a signal with unit power 
the spectral density function for such white noise is 





( _n NG - n - 
^NG < 


(2-16) 


Such a theoretical distribution in single-sided form is also shown in Figure 
2 - 2 . 

It is important to note that the two power spectra shown in Figure 
2-2 are both normalized and thus 


/* DI ( n)dn - / *' DI (n)dn (2-1 7 ) 

-® —CO 

* 1 

The theoretical spectrum Is basically a rectangular pulse function while the 
discrete spectrum is characterized by tapering, The difference between these 
two spectra is generally considered the basis for the second digitization 
effect, 

The preceding descriptions of the two spectra $j(Q) and $ di (q) are 
based purely on mathematical theory. To observe such spectra i n reality the 
corresponding time series would have to be sampled with an infinitesimal 
sampling interval. Actually, finite sampling intervals, T s , must be used but 
this finite (or discrete) sampling process results in aliasing. The aliased 
spectrum, $ + (n), based on the finite sampling process, is related to the 
original spectrum according to the relation [13] 
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(2-18) 


* + (n) = Y *(n+2kn NS ) 

k =-« 

where 

S * Nyquist samp! ing frequency (=»tt/T s ) 


For the case of the discrete white noise 

♦V<» - f 

k=-a» 

x " sin 2 [(JH2kn NS )7r/2n N6 ] 

* — L, 

a NG k=-» t(n+2kn NS )rr/2ft NG J^ 


Numerical evaluation of this series has been carried out for * 100 with 
various ratios of ^ 5 /^, including .5, 1, 2, and 4. The resulting aliased 
spectra are presented in Figure2-3. It is important to note that the figure 
indicates that 


‘•y n> 



(" n NG - a - n NG^ 


(^NG < 


^NG " 


(2-20) 


In this case, by comparision with 

^ l+ Di(^) = ^NG ~ (2-21) 

Thus for white noise the aliasing due to discrete sampling exactly offsets 
the tapering due to discrete generation when the sampling frequency equals 
the generation frequency. Based on this fundamental point, it is clear that 
in the simulation of white noise no tapering of the spectrum occurs as long 
as the sampling rate equals the generation rate. Under most conditions this 
equality is automatically satisfied, 
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The process of convolving the white noise with the appropriate 
impulse response function is also carried out in a discrete manner. The 
process involves selecting (or sampling) values of the white noise signal 
and the impulse response function at equal intervals i n time and then approxi- 
mating the convolution integral by a summation of products. It is important 
to note that the discrete sampling of both the white noise and the impulse 
response function normally occurs at the same rate as the generation rate 
for the white noise. Thus the resulting spectra for the sampled discrete 
white noise, as previously shown, will be uniform, According to the con- 
volution theorem, the spectrum of the output signal equals the spectrum of 
the input white noise multiplied by the product of the Fourier transform of 
the impulse response function and its complex conjugate. Thus, as previously 
noted in subsection 2.4, for a continuous signal, 

* DY (fl) = * DI (n)H(n)H*(n) (2-22) 

The Fourier transform H(ft) for the continuous impulse function, h(t), is 
H(n) = F[h(t)] 

= VT $ DY (n) (2-23) 

The corresponding output spectrum for a discrete signal would be 

*' oy (n) = *' DI (n)H'(n)H'*(n) (2-24) 

The Fourier transform H ' (Q) of the discrete impulse response function, h' (t), 
i s 

H'(a) = F[h'(t)] (2-25) 

Based on the preceding development, for cases in which the sampling 
frequency equals the generation frequency, any difference between the discrete 
turbulence spectrum and the continuous spectrum apparently originates because 
of some difference between H(Q) and H*(&). 
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3. SIMULATED TURBULENCE TAPES 


The turbulence generation procedure described in Section 2 has been 
used to generate six dimensionless simulated turbulence time series which are 
stored cn magnetic tapes as summarized in Table 3-1. The appropriate procedures 
for using the tapes are described elsewhere [1]. Subsection 3.1 provides a 
description of the results of validating the tapes while subsection 3.2 presents 
an explanation of the process for converting from dimensionless to dimensional 
values. 

TABLE 3-1. INDEX OF SHUTTLE 9MLLAIED 
TURBULENCE TAPES (SSTT) 


Tape 

Time Series 

Comments 

SSTT-1 

Uj_ “ gust 

longitudinal gust 

SSTT-2 

U£ " gust 

transverse gust 

SSTT-3 

Uj ■ gust 

vertical gust 

SSTT-4 

- gust gradient 

yaw 

SSIT-5 

9Ug/8Xj " gust gradient 

pitch 

SSTT-6 

auyaxg " gust gradient 

roll 

3.1 

VALDIATION OF SIMULATED TURBULENCE 



A spectral analysis of each of the dimensionless time series has been 
carried out by means of a Fast Fourier Transform FFT4 [ 14 ]. The results, which 
are presented in. Appendix C, demonstrate that the simulated turbulence possesses 
the proper von Karman spectral characteristics. 

All of the dimensionless time series have also been analyzed statis- 
tically to determine the gust and gust gradient probability density functions. 

A shown in Appendix D the results of these analyses indicate that both the 
simulated gusts and gust gradients are normally distributed, with near-zero 
means and standard deviations consistent with the energy content presented 
in Table 2 - 5 . 
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3.2 


CONVERSION TO DIMENSIONAL VALUES 


The dimensionless time series on each tape must be converted to 
dimensional form before actual use in a simulation exercise. The conversion 
process generally involves multiplication and/or division by the appropriate 
turbulence parameters. For dimensionless gusts, u- , the corresponding standard 
deviation, 0^ , should be used. Thus 
* 

u. - o.u. (3-1) 


where 


* 


For dimensionless 
Thus 


dimensional gust 
3u. 

gust gradient, - , the parameters 0. and I. are used, 
dxj n j 



3x 


j 





u.axj 


where 

* 

★ 1 
3u^/3Xj 


* dimensional gust gradient 


(3-2) 


In the case of dimensionless time it is necessary to develop the 
procedures for converting both from dimensionless to dimensional form and also 
to dimensionless from dimensional. In proceeding from dimensionless to dimen- 
sional time the dimensionless time step, T^, represents the basic unit to be 
converted. TTie conversion involves the vehicle velocity, V, and the turbulence 

scale, L.j. Thus 
* 

At. = aL^/V (3-3) 


where * 

At^ * dimensional time step 
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It i s important to note that because both Lj and V vary with altitude, 

the resulting dimensional time step Atj is not a constant. To obtain dimensional 

time, ti, a summation process is involved as follows: 


where 




N 

= rrl At* 
n_1 in 

N 

- al.-E L,„/V. 

Vi " 



L t< Z n> 

V(Z„) 

altitude at nth step 


(3-4) 


In converting to dimensionless fan dimensional time the basic 

unit, the dimensional time step, (St , will normally be a constant. The 

corresponding dimensionless time interval, T. , will be 

* im 


T. = 
im 


v 

aL 


Im 


(3-5) 


The total dimensionless time, t.^, will be 




= El 


m= 1 


im 


M V_5t 

« E m 


m=l aL im 
* M 

- ^-Ev/l 


m*l 


m' im 


(3-6) 
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The dimensionless time, t-^, corresponds to some M' dimension -] ess time inter- 
vals, T.., plus some fractional interval, T\, as follows: 


t 1M - M«T. + V 

where 

0 < T' < T. 


(3-7) 


Thus the number of dimensionless time intervals, M' , can be computed 
as 

M 1 = Int(t. M /T.) 

= Int( lT7 JE V L im ) (3-8) 

where 

Int( ) = integer value of ( ) 

The fractional interval, T', can be computed by the relation 

T = t. M -M«T. (3-9) 

The interpolation process will involve interpolating between t.^, and t^ Vf , ; , 
at the point t-^ as shown in Figure 3-1, 


» — i 

t iM'-l 



Dimensionless 

Time 


Figure 3-1. Relationship Between t-^, t^,, and 


In the conversion to or from dimensional values three parameters 
are required: standard deviation, integral scale o f turbulence, and vehicle 
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speed. The variation of the turbulence standard deviation, a. 
altitude was presented in Table 2-1. The same table contains 
scale, L-, as a function of altitude. The vehicle speed, V, i 
of altitude but also may vary from one trajectory to another, 
provides representative values of V as a function of altitude. 


TABLE 3-2, VARIATION OF SHUTTLE SPEED 
WfTLI ALTITUDE [12] 


ALTITUDE 

(m) 

V 

(m/sec) 

100 

152 

300 

156 

500 

158 

2000 

170 

4000 

188 

6000 

200 

8000 

240 

10000 

300 

20000 

500 

40000 

1928 

60000 

4695 

80000 

1468 

100000 

7521 

120000 

7510 
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4. CIMURM AND RECOMMENDATIONS 


E|- means of a non-recursive discrete generation process, based on 
a vcn Katman spectral model with finite upper limits, dimensionless simulated 
turbulence time series have been developed and stored on six magnetic tapes. 
Longitudinal, transverse, and vertical gusts are simulated as well as the 
gust gradients associated with yaw, pitch, and roll, For each gust or gust 
gradient six separate time series (corresponding to the six altitude bands 
extending from ground level to 120,000 meters) have been stored on each tape. 

The results of spectral analyses of each tape reveals that the 
simulated turbulence possesses the appropriate von Karman spectral character- 
istics, Statistical analyses of the tapes indicate that both the simulated 
gust and gust gradients are normally distributed with near-zero means, Further- 
more the standard deviation of each series is constant with the theoretical 
energy content. 

The Shuttle Simulated Turbulence Tapes (SSTT) are now ready for 
actual use for simulating turbulence at altitudes up to 120,000 meters, 
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APPENDIX A 


DIMENSIONLESS VON KARMAN SPECTRA 
WITH FINITE UPPER LIMITS 


For each altitude band, the three-dimensional spectral model 
for gusts, as given by Eq (2-1), and the three-dimensional model for gust 
gradients, as given by Eq (2-2), have been integrated with respect to fig 
and fij over the finite limits calculated according to Eq (2-3). The six 
resulting one-dimensional spectra are presented in Tables A-l through A-6, 
corresponding to Altitude Bands #1 through 6 respectively, These spectra 
were used in the numerical evaluation of the impulse response functions 
described i n subsection 2.4. 
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TABLE A-l. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 1 
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TABLE A-2. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 2 
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TABLE A-2. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 2 (continued) 
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TABLE A-3. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND I 3 
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TABLE A-3. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 3 (continued) 


CM CO'tr'-fvO'.vfl — r>- 0 ( 7 \ 

cm 'O cm r>- <m tv o ro - t — 

K^M-OON-O-f 
<2 O O - ffl n 00 LO CJ (T, 

r'-«stO’*-roiorM<\itN-* 


05 — '■octocNincM-tfroro 

i7'Min(j)>ONooroM 

•a n «t- in tv in — t in in 
M n 4 N N N N ^ in rf 


— I O 


o » 

I— < 

visa _ 

20 Lil p» 
UJ CO CO 

moj 

M UJ 3 

Q _ I Z 



in cm — o ro ro co ® cd co 

oo-inwvoin®o<i 
o- a-, in in cr. co ro t to <d 
k o ro m m> f- co co co r>. 
IMOOCOiXlOOOOCOOTCDOO 


IN CM CM C\ 
o o o c 
till 
Ul UJ UJ u 

co m ci - cr 

■o CD r- •£ 
n v«io 
•• fS. — 
O'-’On'j 


CM CM <Y CM CM CM 

O o c OOO 

i i i I a t 
wi ui a ui vu ixi 

T ro in ro >t> 

t co n ro cr> o 

eg ro a ro t cj 

ro (v c\ 0 % M> - 9 " 

ro CM CM — 


° >: 

*-• <c 
to 3t oe 

z UJ 
UJ !/) OQ 

zmz 


to O Tf- CT> ro (D CM N •— \D 
co rv in ro eg c? cn r. -d m- 
•=> — cgro-tininvor^co 
ro^ococNinco — m-imo 
ro -o cr- ro -o co ro oo co ro 

oi'tsOiy.-rovoojoro 

— — (M CM CM CM ro ro 




CNiNCMCMcgcM'Nro ro’ro 
oooooooooo 

I I i I i i i I I i 
UJUJUJUIUJUiUJUUJUJ 

CMCDf^^O — CDIOKO 
CM CD 01 M 5 M> O iN tO CM If) 
m-omcctint'Dir 
D'CM — rocD^ — rovoro 
inN-rocg — — — <?■ NM) 


inOMOSWOJCJN - MO 
corMin«cMO(T.r-.xD^- 
o — cm ro m- m in \C* r-- <» 
roMOcncMinoJ — m-n® 
ro vo co ro -43 iti ro -d c- ro 

cnn-sOO' — ro m> © © ro 
_ cm in cm cj ro ro 











TABLE A-4. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 4 














TABLE A-4. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND #4 (continued) 
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TABLE A-5. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND #5 
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TABLE A- 5 . DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND #5 (continued) 
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TABLE A-6. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND #6 
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TABLE A-6. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND #6 (continued) 
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TABLE A-6. DIMENSIONLESS SPECTRUM FOR ALTITUDE BAND # 6 (concluded) 
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APPENDIX B 


ESTABLISHMENT CF LQAER FREQUENCY LIMITS 


The maximum time limit, t max , for which the impulse response 
function is computed, determines the minimum frequency, for which 

the corresponding spectrum is accurately simulated according to the relation 

“ilmin * v/t mx (B-l) 


The simulated turbulence may contain lower frequencies, depending on the total 
length of the time series, but the shape of the actual spectrum for any fre- 
quencies less than will not i n general match the theoretical spectral 

shape. To verify this point for the * 33/11 spectrum two impulse response 
functions were generated, the first extending out to a t fflax of 42.10573 and 
the second extending out to a t max of 85.895689. Two separate turbulence 
time series were then generated, one for each impulse function. The results 
of spectral analysis of these two time series is presented in Figures B-l 
and B-2, As indicated in each figure, the observed spectrum tends to drift 
away from the theoretical spectrum for frequencies below ir/t . The impulse 
response function with the larger t max produces a time series whose spectral 
shape matches the theoretical spectrum to the lower frequency. 


Based on various characteristics of the Space Shuttle simulators, 
a minimum frequency, f lmin* of ,04 hertz was established for generating 
simulated turbulence. The dimensional frequency, ^ i m ^ n * is related to the 
dimensionless frequency, fi -ji m1 - n » according to the relation 




2 * aL i W V 


(B-2) 


Q/ substitution and rearrangement, 

‘max ’ v /(2aL, Amin) 

To satisfy the requirement for a minimum simulation frequency of ,04 hertz, 
values of each impulse response function were computed for 100 dimensionless 
time intervals. As shown in Table B-l, by using a constant number of time 


B-1 


3 





intervals, the actual minimum simulation frequency varies from .0174 hertz 
in Band 1 to .04 hertz in Band 6. Thus, at the lower altitudes the turbulence 
simulation will actually be valid for frequencies somewhat lower than the 
required minimum as given by Eq (B-2). 


TABLE B-l. DIMENSIONAL AND DIMENSIONLESS MINIMUM FREQUENCY LIMITS 





MINIMUM SIMULATION 




FREQUENCY LIMITS 




DIMENSIONAL 

DIMENSIONLESS 

i 

SPECTRA 

BAND 

Mmmomnm 

n ilmin 

1 

*11 


.0174 

.04818 



.0212 

.2176 



| 

.0211 

.3331 




.0396 

.5918 




.04 

7.365 



Mm 

.04 

8.949 

2 

*22’ *22/11 

B 

.0174 

.0206 

.03075 

.2110 




.0211 

.3331 




.0396 

.5918 




.04 

7.365 



HI 

.04 

8.949 

3 

*33’ *33/11 
’*33/22 

B 

.0174 

.0208 

.01865 

.2132 



.0211 

.3331 




.0396 

.5918 




.04 

.4630 



B 

.04 

.5280 


B-4 
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APPENDIX C 


SPECTRAL ANALYSIS CF SIMULATED TUFBULENCE 

B/ means of a Fast Fourier Transform [14] spectral analyses of all 
simulated turbulence have been performed . The results are presented i n 
dimensionless form in Figures C-l through C-36. Table C-l provides a summary 
of these figures. Also included in each figure is the theoretical von Karman 
spectra. The agreement between the theoretical spectra and the computed spectra 
is quite satisfactory. 


TABLE C-l. MATRIX CF SPECTRAL ANALYSIS FIGURES 


SERIES 

TYPE 

ALTITUDE BAND 

1 

2 

3 

4 

5 

6 

U 1 

C-l 

C-2 

C-3 

C-4 

C-5 

C -6 

u 2 

C-7 

C -8 

C-9 

C-10 

C-ll 

C-12 

u 3 

C-13 

C-14 

C-15 

C-16 

C-l 7 

C- 18 


C- 19 

C-20 

C-21 

C-22 

C-23 

C-24 

9Uj/3Xj 

C-25 

C-26 

C-27 

C-28 

C-29 

C-30 

9113 / 8 X 2 

C-31 

C-32 

C-33 

C-34 

C-35 

C-36 


^The spectral analysis involved the first 4096 terms of each time series except 
for bands 5 and 6 for the and Ug gusts. For these cases 8192 terms were used. 


t-1 













a 


Figure C-3. u, - Gust Spectrum, Altitude Band #3 Figubs C 4 - Gust Spectrum, AlPrt^s Etetqd #4 




i 


Figure C-5. Uj - Gust Spectrum, Altitude Band #5 Figure C-6. u, - Gust Spectrum, Altitude 8^d #6 








} 





a 


Figure C-15. u, - Gust Spectrum, Altitude, Band #3 Figure C-16 u 3 - Gust Spectrum, Altitude Band #4 



Figure C-17 ^3 -oust Spectrum, Altitude Band #5 Figure C-18. u- - Gust Spectrum, Altitude Band #6 




c-u 


Figure 019. 3U2/3XJ - Gust Gradient Spectrum, Figure 020. aug/SXj - Gust Gradient Spectrum, 

Altitude Band #1 Altitude Band #2 










a 


Figure C-25 3 uj/3Xj - Gust Gradient Spectrum, Figure C-26. au^/ax^ - Gust Gradient Spectrum, 

Altitude Band #1 Altitude Band #2 











a 


fgure C-33. 31 * 3 / 8 X 2 - Gust Gradient Spectrum, Figure C-34. au^/axg - Gust Gradient Spectrum, 

Altitude Band #3 Altitude Band #4 



C-19 


Figure C-35. 3u 3 /3x 2 - Gust Gradient Spectrum, Figu^ C-36. au 3 /gx 2 - Gust Gradient Spectrum, 

Altitude Band #5 Altitude Band #6 


APPENDIX D 


STATISTICAL ANALYSE OF SIMUL ATED TURBULENCE 

Ey means of standard statistical analysis procedures each of the 
SSTT has been analyzed to determine its mean value, standard deviation, and 
probability density distribution . The resulting mean values are presented 
in Table D- 1 while Table D-2 contains the resulting standard deviations. 

A expected all mean values were near zero. The standard deviations represent 
the square root of the energy content. The ratio of the theoretical energy 
content (from Table 2-5) to the square of the corresponding standard deviation 
(from Table D-2) is presented in Table D-3. The agreement appears quite satis- 
factory . 

The gust and gust gradient probability density distributions are 
presented in Figures D-l through D-36 in accordance with Table D^l. In each 
figure the corresponding theoretical normal distribution is also presented. 

The results indicate that both the gust and gust gradient time series are very 
close to normal distributions. 

TABLE D-l. MEAN VALUE OF GUST AND GUST CRACJENIS 


SERIES 

ALTITUDE BAND 

TYPE 

1 

2 

3 

4 

5 

6 

U 1 

-.019295 

-.042050 

-.051852 

-.088142 

-.0455 

-.0464 

u 2 

-.010671 

-.029576 

-.0371 

-.049431 

-.0428 

-.0441 

u 3 

-.006806 

-.029652 

-.037043 

-.049370 

-.043788 

-.046637 

3U2/8XJ 

-.000002 

-.001572 

-.002794 

-.005628 

-. 172152 

-.206385 

8113/ ^2 

-.000001 

-.001591 

-.002798 

-.005649 

-.004293 

- ,004893 


-.005760 

-.073072 

-.103823 

-,160303 

-. 171448 

-.288178 


The statistical analysis involved the first 4096 terms of each time series except 
for bands 5 and 6 for the and U£ gusts. For these cases 8192 terms were used. 

D-1 


1 




TABLE D-2. STANDARD DEVIATION OF GUST AND GUST GRADIENTS 


SERIES 

ALTTTUDE BAND 

TYPE 

1 

2 

3 

4 

S 

6 

U 1 

.788959 

.927098 

.946351 

.96427 1 

.99888 

.99996 

u 2 

.707845 

.925201 

,94619 

.964152! 

.99764 

.99863 

u 3 

.524571. 

,915606 

,938552! 

-958985 

.961845 

.967651, 

3u2/3Xj 

.76662 7 

3.625937. 

4.976618' 

7.3568081 

41.717292 

48.0527341 

3u 3 /3x2 

.390512 

3.488677 

4.758426 

7.037516 

6.458092 

7.2166481 

3u 3 /3Xj 

.394539 

3.5082801 

4.784378 

7.075349 

9.778736 

19.790119) 


TABLE D-3. RATIO OF THE THEORETICAL ENERGY CONTENT 

TO THE SQUARE OF THE OBSERVED STANDARD DEVIATION* 


EH 

ALTTTUDE BAND 

ill 

1 

2 

3 

4 

5 

6 

U 1 

1.0001 

1.0000 

1.0000 

1.0000 

.9999 

1.0001 

u 2 

.9999 

1.0000 

.9999 

1.0000 

1.0000 

1 0000 

u 3 

1.0001 

1.0000 

1.0000 

1.0001 

1.0000 

.9999 

3u2/3Xj 

1.0000 

1.0000 . 

1.0000 

1.0000 

.9998 

1.0000 

3u 3 /3x 2 

1.0000 

1.0000 

1.0000 

1.0000 

1.0001 

1.0000 

3u 3 /3Xj 

1.0002 

1.0000 

1*0000 

1.0000 

1.0000 

.9999 




Theoretical energy content taken from Table 2-3. 

'Observed standard deviation taken from Table D-2. 
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TABLE D-4. MATRIX OF STATISTICAL ANALYSIS FIGURES 


SERIES 

TYPE 

ALTITUDE BAND 

1 

2 

3 

4 

5 

6 

U 1 

D-l 

D-2 

D-3 

D-4 

D-5 

D-6 

u 2 

D-7 

D-8 

D-9 

D-10 

D-ll 

D-12 

u 3 

D-13 

D-14 

D- 15 

D-16 

D-17 

D - 1 8 

ai^/SXj 

D- 19 

D-20 

D-21 

D-22 

D-23 

D- 24 

3u 3 /3Xj 

D-25 

D-26 

D-27 

D-28 

D-29 

D-30 

9u 3 /3x2 

D-31 

D-32 

D-33 

D-34 

D-35 

0-36 
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Figure D~7 u 2 - Gust Probability Density Distribution, Altitude Band #1 
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D-12 


Figure D-9. - Gust Probability Density Distribution, Altitude Band #3 
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Figure D-16. u^, - Gust Probability Density Distribution, Altitude Band #4 



o 



Q 


-Q 

O 

i- 

<x 


CO 

13 

CD 


CO 

3 


I 

o 

<u 

s- 

13 

o> 




D-20 


DIMENSIONLESS OUST STRENGTH 



A1ISN3Q AimayaoHd 

D-21 


i 


DIMENSIONLESS G-UST STRENG-TH 


o 

=> 



D-22 


4 


Figure D-19. 3u 2 /3xj - Gust Gradient Probability Density Distribution, Altitude Band #1 
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Figure D~22. 3u ? /3x, - Gust Gradient Probability Density Distribution, Altitude Band #4 
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Figure D-24. SUg/SXj - Gust Gradient Probability Density Distribution, Altitude Band #6 
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Figure D-25. 3u.,/3x, - Gust Gwudient Probability Density Distribution, Altitude Band #1 
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Sfgure D-29 3ia 3 /3xj - Gust G^disnn Probability Density Distribution, Altitude Bu,nd #5 
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Figure D-32. 3u^/3x2 - Gust Gradient Probability Density Distribution, Altitude Band #2 
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